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ABSTRACT 

At the epoch of reionization, when the high-rcdshift inter-galactic medium (IGM) is being enriched 
with metals, the 63.2/im fine structure line of 01 is pumped by the ~ 1300A soft UV background 
and introduces a spectral distortion in the Cosmic Microwave Background (CMB). Here we use a toy 
model for the spatial distribution of neutral oxygen, assuming metal bubbles surround dark matter 
halos, and compute the fluctuations of this distortion, and the angular power spectrum it imprints on 
the CMB. We discuss the dependence of the power spectrum on the velocity of the winds polluting 
the IGM with metals, the minimum mass of the halos producing these winds, and on the cosmic epoch 
when the OI pumping occurs. We find that, although the clustering signal of the CMB distortion 
is weak ((Sy) rms < 1CP 7 ) (roughly corresponding to a temperature anisotropy of ~ nK), it may be 
reachable in deep integrations with high-sensitivity infrared detectors. Even without a detection, these 
instruments should be able to useful constraints on the heavy element enrichment history of the IGM. 
Subject headings: cosmology: cosmic microwave background - cosmology : theory - galaxies: inter- 
galactic medium - atomic processes 



1. INTRODUCTION 

Probing the Dark Ages - the epoch between the last 
scattering surface of the CMB at z ~ 1000 and the com- 
pletion of the reionization of the IGM at z ~ 6, including 
the formation of the first luminous objects - constitutes 
the next frontier of observational cosmology. Among the 
major open questions are the nature of the objects that 
rcionized the universe, and the origin of the first heavy 
elements, as well as the efficiency with which they were 
mixed into the high-redshift IGM. 

To date, the reionization history has been con- 
strained by observations of Ly man series abso rption 
spectra of z ~ 6 quasars (see iFan et akl [2006 for a 
recent review ) by GRB observations at high redshifts 
(jTotani et al.l I2006D . by observations of Lyman-q se- 
lected galaxies (iKashikawa et al.l 120061: iMcQuinn et al.l 
l2007HDiikstra et al.ll2007tlMalhotra fc Rhoadsll20q6D and 
by CMB polarization anisotropics (|Kogut et alj 120031 : 
iPage et aTTl2007HSpergel et al.l l2007). The present con- 
straints, however, are consistent with a wide range of 
scenarios, and considerable theoretical and experimental 
effort is devoted to developing new probes of the high- 
redshift universe. A promising possibili ty is to use the 
redshifted 21cm hyperfine line of HI (see lFurlanetto et al.l 
( 2006) and references therein) . Metal enrichment of the 
high-redshift IGM may be anot her useful tracer o f the 
reionization process. For example | Basu et alj (|2004f ) and 
iHernandez-Monteagudo et a l. ( 2006) considered the elas- 
tic resonant scattering of CMB photons by intergalac- 
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tic metals. Recent studies have focused, in part i cular , 
on the detectab ility of neutral oxygen (OI; lOhl 120021 : 
iBasu et aT1l2004l ). as this element is t hought to be pro- 
duced in abundance by the first stars (jHeger fe Wooslevl 
2002). Scattering of UV photons by OI, and the corre- 
sponding absorption features in the spectra of quasars - 
the OI for est- was proposed as a possible observable by 
loU(l2002h. and may have recently been detected at z ~ 6 
(|Becker et aT7~ll2006f) . Since OI and HI are in charge ex- 
change equilibrium, oxygen is likely to be highly ion- 
ized in regions where hydrogen is actively ionized. How- 
ever, the recombination time for oxygen is shorter than 
the Hubble time, and it can be neutral even in "fossil" 
HII regions where hydrogen had been ionized, but where 
short-lived ionizing s ources hay e turned off, allow i ng the 
region to recombine (fOhll2002f) . I Oh fe Haimanl (|2003l ) 
show that the filling factor of such fossil HII regions can 

be large ( <; 50%) prior to reionization. 

In a previous paper (|Hernandez-Monteagudo et al.l 
l2007t hereafter Paper I), we showed that the 63.2^tm 
fine structure line of neutral OI can be pumped by the 
- 1300 A soft UV background, via the Balmer a line 
of OI. This is analogous to the Wouthuysen-Field ef- 
fect for exciting the 21cm line of cosmic HI. In Paper I, 
we found that OI at redshift z should be seen in emis- 
sion at (1 + z)63.2/im, and for 7 < z < 10 it would 
produce a mean spectral distortion of the CMB with 
a y-parameter of y = AI„/B u [Tcmb] = (10~ 9 — 3 x 
10- 8 )(Z/10- 3 Z Q )(4rv7lO), where Z is the mean metal- 
licity of the IGM and Ijjv is the UV background inten- 
sity at 1300A in units of 10~ 21 erg/s/Hz/cm 2 /sr. In 
principle, this distortion could be detectable through a 
precise future measurem ent of the CMB spectrum (see 
iFixsen fc Matherl |2002 for the prospects of such mea- 
surements), and would then open the possibility of per- 
forming tomography of the metal distribution. In com- 
bination with HI 21cm studies, it could yield direct mea- 
surements of the abundance and spatial distribution of 
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metals in the high-redshift IGM. 

Since oxygen pollution at high-rcdshift is likely associ- 
ated with most over-dense regions in the IGM, hosting 
the first star forming activity, there must be spatial fluc- 
tuations in the 01 abundance, causing fluctuations in the 
corresponding y distortion. In this paper, we consider 
the clustering properties of this signal. In particular, we 
use toy-models to describe the metal distribution, and 
compute the two-dimensional angular power spectrum 
of the CMB intensity. This would be appropriate for an 
instrument with a single (or a few discrete) frequency 
bands that probes the metals in a single (or a few dis- 
crete) narrow redshift bins. We find that the cluster- 
ing signal induced by the inhomogeneous 01 pumping is 
small, but possibly detectable with forthcoming instru- 
ments, e.g., with deep integrations (for 2-3 months) with 
the full ALMA array. A detection of the metal enrich- 
ment during the Dark Ages would be complementary to 
HI 21cm studies, and would provide clues about the dis- 
tribution of metals in the IGM before galaxy formation 
started in full. 

The rest of this paper is organized as follows. In §[2j wc 
discuss the basics of the 01 pumping process. In § [3l we 
outline the computation of the distortion induced by 01 
pumping along a given line of sight in the CMB. In §[4] we 
discuss our method to compute the angular power spec- 
trum of the distortion, with the basic assumption that 
metals trace the distribution of collapsed halos. In § [5l 
we introduce toy models to describe the spatial distri- 
bution of metals, with the basic assumption that metals 
cluster around dark matter halos. In § [5J we present 
our main result on the angular power spectrum, and dis- 
cuss its dependence on the basic model parameters, as 
well as its dctcctability. Finally, in §0 we summarize 
the implications of this work and offer our conclusions. 
Throughout this paper, we adopt a set of "concordance" 
cosmological parameters for a flat universe, 
SIa = 0.71, Slf, = 0.047, h = 0.72, with a power spectrum 
normalization a% = 0.75 and slope n — 0.99. 

2. THE BALMER-a PUMPING OF OI 

In Paper I, we computed the distortion in the CMB 
induced by neutral oxygen in environments where UV 
radiation is pumping the fine structure 63.2/im Ml tran- 
sition (between the n = 2 electronic states 3 P2 and 3 Pi) 
via the Balmer-a line at ~ 1302A connecting these two 
states with the excited n = 3 electronic state 3 Si. In 
what follows, we shall denote the 3 P2 and 3 Pi states 
as "0" and "1" respectively, and the excited state 3 Si 
as " 2" . Energy and frequency differences between levels 
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will be denoted as E 



, with i,j = 0, 1, 2 and j > i. 
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In Paper I, we showed that this process of Balmer a 
pumping modifies the occupation of the levels and 1 
and therefore introduces a shift in the spin temperature 
Ts, so that it slightly departs from the CMB temper- 
ature, Tcmb- Since at reionization the UV background 
flux at viq is smaller than at v%\, Ts will be above Tcmb, 
producing an excess of 63.2 /im photons. Each of these 
excess photons is generated by the sequence of three ra- 
diative transitions — > 2 — * 1 — > 0, i.e., a photon of 
frequency ^20 is broken into two photons of frequencies 
v 2 i and v w . 

In Paper I, we were interested in the average distortion 
over the full sky, and considered the optically thin limit 
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Fig. 1. — The distortion parameter y induced by the Balmer-a 
pumping of the OI 63.2 (im transition is shown as a function of the 
optical depths no and T20 at redshift 2 = 7. The optically thin 
solution presented in Paper I is shown as the solid curve, accurately 
describing the system at low and intermediate no- When the 2<->0 
transition becomes optically thick (T20 ~ 1, corresponding to no ~ 
3 X 10~ 5 ), the results are obtained by tracking the photon and level 
populations as dictated by equation and shown by the dashed 
curve. Note that at low no there is a slight mismatch between the 
solid and the dashed lines, caused by numerical error due to the 
finite size of our integration step. 

for the Balmer-a transition. While this is justified if the 
OI is homogeneously distributed (in which case, T20 ~ 
O.O1[Z/(1O~ 25 Z )] and r 2 i < 1), here we will consider 
metal distributions that can be highly inhomogeneous. 
The metal-enriched bubbles can potentially contain a 
large concentration of metals and become optically thick 
at 1^20 ■ For this reason, we consider the case T20 <; 1 and 
self-consistently solve the equations of radiative transfer 
to follow the UV intensity as the background radiation 
penetrates the OI-rich bubble. 

The equations governing the evolution of the relative 
populations rij (with i = 0, 1) and the number of photons 
at the frequencies riff , j = 1, 2 are: 



dno 
~dt 



-no 
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+ P01 
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dnf = Xl 
dt 8tt l 



~P in + P w ni 



(1) 
(2) 



where ^ v denotes the line profile, assumed to have a 
width corresponding to a thermally broadened line with 
Toi ~ 20 eV (corresponding to a relative thermal width 
of ~ 4 x 10~ 5 ). For simplicity, we assume top-hat shapes 
for the line profiles, and we include the effect of cosmolog- 
ical redshift by adopting an appropriate effective length 
for the line-of-sight integral (discussed below). In addi- 
tion, atom and photon number conservation imply that 
dni/dt = —dno/dt and dn 2 J- jdt — —dnfP/dt. 
As initial conditions, 



2 x 10 20 erg cm 2 s 



we chose T s = Tcmb-, tf 1 
11/ 1 sr" 1 , and i 20 // 21 
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with j3 = 0.02 as in Paper I (our results are largely in- 
sensitive to the choice of this parameter, as shown in 
Paper I) . Figure [1] shows the result of evolving this sys- 
tem of equations at z = 7. At low t\q, the results (dashed 
curve) coincide with the optically thin steady state solu- 
tion of Paper I (solid curve) within a few percent level 
corresponding to innacuracies in the numerical integra- 
tion. At larger values of t\q <; 10 -4 , the Balmcr-a transi- 
tion becomes optically thick and the number of available 
pumping UV photons decreases. As a result, once these 
photons are exhausted, the distortion y flattens out at 
the maximum value of y ~ 3 x 10 -7 . 

3. THE SPECTRAL DISTORTIONS ON THE CMB 

The total distortion introduced by 01 in the CMB at 
a given observing frequency v obs is 



V= / dr(j v /B v [T ])(r) 



3v 



(v ohs )Ar, (3) 



where j„ is an effective emissivity which is non zero only 
at the redshift satisfying 1 + z s « v\ol v obs (see Paper I 
for more details). For a uniform 01 distribution, and as- 
suming that the 01 atoms follow the Hubble expansion, 
the effective length Ar along which this emissivity is not 
zero is cH' 1 (z s ) (A^/j/^io), where Ai/th is the width 
of the line. However, if oxygen is clumped (for exam- 
ple, confined within bubbles) then the effective length 
Ar may vary. If the physical size L of the bubble is 
small, L <> cH^ 1 (z s ) (Ai/ t ^/i/ 10 ), then Ar « L. On the 
other hand, if L ^> cH^ 1 (z s ) (Av t h/vio) and the bubble 
has split from the Hubble flow, then also Ar = L. But 
if the bubble is comoving with the Hubble flow, then 
Ar = cH~ 1 (z s ) (Avth/vio), (this length, at z = 7, corre- 
sponds to ~ 12 Kpc). Note that, unless otherwise stated, 
distances will always be in physical units (and not co- 
moving). In this work, we shall assume that bubbles are 
not gravitationally bound, and for simplicity, we further 
assume that they expand with the Hubble flow, so that 
Ar = min (L, cH~ 1 (z s ) Aiy t h/^io)- Hence, in this case, 



y ■ 



B U [T 



min I L, cH 1 (z s ) 



(4) 



Depending on the angular and spectral resolution of a 
specific instrument, the y-distortion along a given line- 
of-sight, and at a given frequency, may have to be con- 
volved with the respective response functions (clearly, 
such smoothing can reduce the fluctuations). In what 
follows, we adopt a model instrumental beam (or PSF) 
response function V^ s /(fi) and a frequency response 
function $(v). Both will be taken to have accurately 
determined Gaussian shapes, so that J dhV ps f{h) = 
J dv<fr(v) = 1, where n denotes the unit vector. 

We assume generically that each OI-bubble is associ- 
ated with a halo of some mass M, with an abundance 
given by the Sheth-Tormen (ST, ISheth fc Tormenll2004h 
mass function dn/dM, and that the 01 atoms are dis- 
tributed in a bubble around the halo with a spherically 
symmetric density profile W b (r); in this second 
convolution, over the bubble density profile, is necessary. 

If we neglect the peculiar velocities and internal mo- 
tions of the bubbles, the effective distortion along the 



line of sight in direction n, and at frequency v bs, reads 

roc /> / i r s +Ar/2 

y e //Kbs,n) = / / dvdxi$(v)V ps f(ii) / dr 

JO J PSF Jr s -Ar/2 

^M^(y,M)^(y-r)^ (5) 

where r = rn. The distance r is centered at the 
resonant value r s = r[z s ] and Ar corresponds to the 
effective length defined above. The integral over M 
counts 01 bubbles around halos of various masses (above 
some minimum halo mass M min , as discussed further 
in § [5] below), whereas the integrals over v and n de- 
scribe the cosmological volume probed by the PSF and 
the spectral response of the detector. Let us change 
variables from frequencies and angles to spatial coor- 
dinates, and define the three-dimensional instrumental 
PSF, B(x-r.) = <b{v[x]){v obs /rl){dz/dr{z))V psf (nf. 
Further defining y(M,z) = (j v /B v [T ])Ar, the distor- 
tion generated by a bubble around a single halo of mass 
M at redshift z, we can re-write equation ([5]) more trans- 
parently as 



y e //Kb S ,n) 



dM 



( dn 
B *{dM* W » 



V- (6) 



where * denotes convolution. In what follows, the twice 
convolved halo number density will be denoted by n, i.e., 
n(M, r s ) = [B-k (dn/dM * W b )} (M, r s ) 

Finally, the physical collision between two metal-rich 
bubbles will lead to a complicated thermo-dynamical in- 
teraction. Rather than modeling this process, the above 
prescription assumes a linear addition of the 01 density 
(or y-distortion) from two physically overlapping metal 
bubbles. This simple assumption at least captures the 
enhancement of the metallicity in the overlap regions. 
We note that this is different from the case of merging 
HII bubbles during reionization, when mergers conserve 
volume and result in the expansion of the joint bubble. 

4. THE EFFECT OF CLUSTERING 

We next compute the second order moments (corre- 
lation function and power spectrum) of the distortion 
field generated by 01 at 63.2/^m during reionization. 
Throughout this discussion, we will consider a fixed red- 
shift, and suppress the frequency-dependence of y e //(n). 
The angular correlation function can be written as 



(y e //(niW/(n 2 )) 



dM 1 dM 2 (nin 2 )yiy2 



= I dM 1 dM 2 nxhiVih ( l + &h(ri,Mi,r 2 ,M 2 ) 

V) 

In this equation, h denotes average n. The correla- 
tion function £hh{ r i, Mi, r 2 , M 2 ) corresponds to the halo- 
halo correlation function convolved with both the window 
function of the profile of the 01 distribution (Wb) and the 
window function of the experiment (B). Keeping in mind 
that the Fourier counterpart of ^hh is the matter power 

5 The integral f dxB(x — r s ) = 1 is normalized to unity in the 
volume centered at r s . 
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spectru m times the square of the bias factor b(M, z) com- 
puted in lSheth fc Tormeil (|1999f ). equation ([J|) becomes 



eff 



(ni)y e //(n 2 )) = / dM x dM 2 fiifi 2 yiy 2 hb 2 X 



dk 



w, 



b.k 



exp(-ik(ri - r 2 )), 



with P m (k,z s ) the matter power spectrum at redshift 
z s , (note that we have dropped a constant -k = 0- term 
and hence we are looking at the departure of the corre- 
lation function from its mean value) and W b ^ and B-^ 
the Fourier counterparts of Wb and B respectively. Since 
we are assuming that every halo is producing a bubble, 
and the dominant signal will be due to the clustering 
of different bubbles, we can neglect non-linear correc- 
tions to the matter power spectrum. Such non-linear 
corrections would boost the clustering signals we predict 
below, but o nly on comoving scales of ^0.1 Mpc at 
z « 6 (e.g. Illiev et al. 2003), which is well below the 
scale at which the 01 clustering signal peaks. The inte- 
gral on k can be split into a transverse (kj_) and a par- 
allel (k z ) component along the line of sight. We assume 
that both B^, W b ^ can be factorised as B^ — (B^ Bu z ) 
and W b k = (W b \ l± Wb,k z )- For a Gaussian instrumental 
response, we can write = exp — [kj_(r s ab) 2 /2] and 

Bu z = cx-p — {kl((dr(z)/dz)a u /v bs) 2 /2}, with o b related 
to the width of the angular PSF and cr„ to the width of 
§(y). For simplicity, the bubble is taken to have a Gaus- 
sian profile and a volume equal to (\/2ttL) 3 . After this 
decomposition, the correlation function reads 



dkj_ 



(y e //(ni)2/ e //(n 2 )) : 
exp (-ikx(ri - r 2 )j_) / dM 1 dM 2 nih 2 yiy 2 hb 2 



dk 

\B k ,W b ,k, | P m (k±, k z ,z rs ). (9) 

Z7T 



Finally, by noting that in the flat sky approximation 
k^ ~ we obtain the angular power spectrum: 



dM x dM 2 hih 2 yiy 2 b x b 2 

J \^W b , k fP m {\/r s ,k z ,z s ). (10) 

Let us briefly examine the behaviour of the C/'s. We 
first recall that the Fourier window function of a bub- 
ble is proportional to the bubble volume (oc L 3 ), and 
that the oxygen_ number density in bubbles is propor- 
tional to no 1 1 (NL 3 ), with N the average bubble num- 
ber density, noi the global mean oxygen number density 
at a given redshift, and L the typical bubble size. We 
also note that the bias factor b almost cancels the red- 
shift dependence of the growth factor of perturbations 
(e.g.. lOhTe t al. 2003j), so the scaling of the maximum of 
l 2 Ci is {l 2 Ci) max cx Q 2 n 2 OI (Ar) 2 , with Q a frequency- 
dependent function accounting for the efficiency of the 
pumping process and Ar the effective length in equa- 
tion ([3]). Q dominates the overall redshift dependence 
of the power spectrum amplitude. If L is smaller than 
cH- 1 (z s ) Avth/vw, then (l 2 Ci) max oc Q 2 n 2 OI i 2 , but 
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Fig. 2. — The radii of metal bubbles around halos of mass M 
produced in our fiducial toy-model, assuming metals expand over 
the typical age of the halo at a constant speed of (v^ — »esc) 2 - 
Here vi, is a wind velocity equal to 150km/s, and the escape 

velocity from the halo. The solid curve shows the radii at z = 6.5; 
the dashed curve at z = 7.5. The sharp feature at M ft: 10 9 Mq 
corresponds to the limit where metals are trapped inside the halo 
potential well (vi, as v esc ). The bubble radius for more massive 
halos is assumed to be a fixed fraction (here taken to be 10%) of 
the virial radius, and so it scales as <x M 1 / 3 . 

otherwise (l 2 Ci) m ax oc Q 2 n 2 Q1 , i.e., the band power 
spectrum is independent of the bubble size (except that 
the bubble Fourier window function suppresses power 
at scales smaller than the bubble size). In either case, 
note that (l 2 C'i) max is proportional to the square of 
the oxygen metallicity. However, if the OI abundance 
within the bubbles is so high that they become optically 
thick at J/20, then_y will reach its plateau, y ~ Q and 
{l 2 Ci) max oc Q 2 (NL 2 ) 2 , i.e., {l 2 Ci) max will be propor- 
tional to the square of the number density of bubbles 
times their cross-sectional area, assuming that the fre- 
quency resolution of the instrument corresponds to a ra- 
dial distance that does not exceed the typical bubble size. 
Otherwise the scaling would be (l 2 Ci) max oc Q 2 (NL 3 ) 2 . 
The discrete nature of the source distribution means 
that in addition to the correlation term computed above, 
there will be a Poisson term. The Poisson contribution 
can be computed by considering the limiting case of equa- 
tion (|10p for a random source distribution and an exper- 
iment of infinite angular and spectral resolution, 

Ci ~ y 2 {NL 3 )AQ b (11) 

with y the average bubble distortion, NL 3 the volume 
fraction occupied by bubbles and AQb their typical an- 
gular size. In practice, we find that the shot noise con- 
tribution is nearly always sub-dominant in the results we 
present below. 

5. TOY MODEL FOR METAL DISTRIBUTION 

To explore the detectability of the angular power spec- 
trum of the OI-induced CMB distortion, we need a 
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model for the spatial distribution of metals in the IGM. 
Following the discussion in §[3] and §|4] above, this means 
we need to specify the profile W&(r) around a halo of 
mass M at each redshift z. 

Our basic simplifying assumption is that metals are 
confined in bubbles that had previously expanded into 
the IGM, but that these bubbles have settled to follow 
the Hubble flow. We assume that when a halo is formed, 
a metal-polluting wind is launched from the galaxy at 
the center of the halo, producing a bubble that subse- 
quently expands at a constant velocity for the age of the 
halo. In principle, halos of the same mass that exist at 
a fixed redshift can have a distribution of ages. For sim- 
plicity, however, we assign a fixed age to halos of mass 
M at redshift z, corresponding to the average redshift Zf 
at which such halos assembled 65% of th eir mass (see, 
e.g., equation 2.26 in lLacev k. Cold ri993f l. At redshift 
z, the radius of the bubble around a halo of mass M is 
therefore given by 

L(M, z) = 0.1rw(M, z f ) + (v 2 - v 2 esc ) x l 2 [t{z) - t(z f )] , 

with r v i r the virial radius of the parent halo, and t(z) 
the cosmic time at redshift z. The first term on the right 
hand size represents a rough estimate for the size of the 
galaxy in the halo (i.e., the region producing the met- 
als). The velocity v b , taken to be a constant in the range 
of 50 — 1500 km s , represents the wind velocity "in 
vacuum". We then subtract u esc , corresponding to the 
escape velocity from the parent halo at a distance of 0.1 
r V i r (M, Zf), from this wind velocity, to take into account 
the fact that metals loose energy as they travel out of the 
gravitational potential well. Since we assume that met- 
als have no peculiar velocities, the bubbles start comov- 
ing with the Hubble flow by the time they are observed. 
Equation (|12[) assumes further that star-formation takes 
place instantly when the halo forms. In Figure[2j we show 
the bubble sizes predicted under these assumptions. 

For simplicity, we normalize the total amount of oxy- 
gen mass in all bubbles at a given redshift by specifying 
the global average metallicity of the IGM at each red- 
shift, i.e., 

Z av (z)= dM-£-(M,z)Z(M,z)V b (M,z), (13) 

JM mln dM 

with Z(M, z) the metallicity within a bubble at red- 
shift z originated in a halo of mass M and with volume 
Vb = 4ir/3L 3 , and dn/dM is the halo mass function. For 
halos with Vb > v esc , we further assume that, at a given 
redshift, the oxygen abundance within the bubble is pro- 
portional to the mass of the parent halo. Provided that 
the bubble radii are practically constant in this mass 
range, this is roughly equivalent to assuming that the 
star formation rate scales linearly with halo mass, and is 
constant throughout the age of the halo. For halos with 
Vb < v esc , we assume that the oxygen density is constant 
(since the metals produced in these halos are assumed to 
be trapped in a volume Vb oc M); however, these larger 
halos contribute negligibly to the clustering signal below 
and our results are insensitive to this assumption. 

In Figure [3j we show the evolution of the global av- 
erage (dashed lines), and bubble- volume weighted (solid 
curves) metallicities versus redshift. The dotted curves 
show the fraction of the total volume contained inside 



the bubbles (note that metal-rich bubbles can overlap 
and this fraction can exceed unity). In all panels, the 
global average oxygen abundance is scaled to 10~ 2 ' 5 Z Q 
at each redshift. This value is not unrealistic, as there 
are observational constraints on the cosm ic metallicity 
indic ating that Z/Z Q ~ 10 -3 at z > 4 (|Schave et alJ 
2003) , and oxygen from first-generation stars is likely to 
be somewhat overa bundant relative to the solar value 
(jMevnet et all l2006h . The highest redshift, z = 12, 
shown in the figures represents the earliest epoch when 
this floor metallicity may have been established (e.g. 
lHaiman fc Loeblll997f l. Note that since the global metal- 
licity is re-scaled at each redshift by hand, the curves in 
Figure should be interpreted as independent models 
of the metal distribution at each redshift, and not as 
evolutionary models. In order to illustrate the metal dis- 
tribution under different assumptions, in panel (a), we 
assume M min = 5 x 10 5 M Q and v b = 50 km/s; panel (b) 
is the same as (a), except with Vb = 150 km/s, while in 
panel (c), we assume v\, — 150 km/s and M m j„ = 5 x 10 7 
M Q . 

The enrichment scenarios considered in this work are 
comparab le with those prese n ted in previous studies. In- 
deed, in iFurlanetto h Loebl (|2003f) the typical bubbles 
sizes at z ~ 6 are a few tens of kpc, in agreement with our 
estimates, whereas, after looking at the simulations of 
lOppenheimer etltfl (|2007h , we have found similar enrich- 
ment le vels in the relevant redsh ift ranges (z <G [5,10]). 
Only in lScannapieco et al.l (|2002l ) the bubble filling fac- 
tors tend to be slightly smaller than in our case, since, 
motivated by WMAP observations, we are considering 
early reionization models. 

6. RESULTS AND DISCUSSION 

The angular power spectrum (correlation term plus 
shot noise term) of the CMB distortion produced by the 
OI pumping is shown in Figure [4] We have adopted 
an angular resolution of one arc-second, and a relative 
spectral resolution of Av/v — 10 -4 . The angular reso- 
lution corresponds to I ~ 7 x 10 5 and therefore should 
resolve all bubbles under consideration. The frequency 
resolution corresponds at z = 6.5 to roughly 60 kpc in 
physical units. In all three panels, we consider varia- 
tions from our adopted fiducial model given by z = 6.5 
(observable at 630 GHz), I 20 = 6 x 10- 20 erg cm" 2 s" 1 
Hz" 1 sr" 1 , M min = 5 x 10 5 M Q , and v b = 150 km/s, 
corresponding to L ~ 60 kpc at that redshift. Panel 
(a) shows the dependence of the signal on bubble expan- 
sion velocity/bubble size: 150 (solid curve), 60 (dashed 
curve) , and 200 (dotted curve) km/s, which correspond 
to I ~ 60, 23 and 80 kpc respectively. If Vb = 60 
km/s (L ~ 23 kpc), then L < cH~ x (z s )Av t h/ '^io, and 
the bubbles remain optically thin (y has not reached its 
maximum value), so that (l 2 Ci) max cx Q 2 L 2 . On the 
other hand, for U(,^150 km/s (L ~60 kpc), we are in the 
regime of L >cH^ 1 (z s )Av th /v w and {l 2 Ci) max cx Q 2 . 
Therefore, if Vb = 200 km/s, the amplitude of the spec- 
trum should not change noticeably. However, the bubble 
size (L ~ 80 kpc) is slightly larger than in our fidu- 
cial model (L ~ 60 kpc), and when comparing both 
cases one can clearly see the effect of the bubble win- 
dow function suppressing power at small scales. There- 
fore {l 2 Ci) max oc (i>{,/150 km/s) 2 if Vb < 150 km/s. In 
panel (b), we find a strong dependence of the power on 
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Fig. 3. — Models of metal enrichment in the IGM at various redshifts. The global average metallicity Z/Zq is shown by the dashed 
curves, and the metallicity inside individual metal bubbles by the solid curves. Note that the global metallicity is renormalized at each 
redshift by hand to the same metallicity, and the curves should not be interpreted as evolutionary models. The dotted curves show the 
volume-filling fraction occupied by bubbles (note that metal-rich bubbles can overlap and this fraction can exceed unity). Panel (a): 
Mmin = 5 X 10 5 Mq, Vf, = 60 km/s and global average Z/Zq = 10 — 2 5 set to a constant in redshift. At z ~ 6.5 the typical bubble size is 
L ~ 22 kpc, and bubbles have been generated in halos of typical mass M ~ 4 X 10 6 M . Panel (b): same as panel (a) but = 150 km/s, 
which corresponds, at z ~ 6.5, to a typical bubble size of L ~ 60 kpc. Panel (c) vt, = 150 km/s (corresponding to L ~ 50 kpc at z ~ 6.5) 
and M min = 5 x 10 7 Mq. 
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Fig. 4. — The dependence of the angular fluctuation power, l 2 Ci, on metal distribution model parameters and on redshift. In all 
three panels, we display our adopted fiducial model given by: z = 6.5 (observable at 630 GHz), I 2 " = 6 X 10 _20 erg cm -2 s _1 , Hz -1 
sr _1 , M m i n = 5 X 1O 5 M0 and L ~ 60 kpc. Panel (a): dependence of the signal on bubble size: 23,60,80 kpc for the dashed, solid and 
dotted curve, respectively. For L ~ 23 kpc we still have L < cH~ 1 (z a )Au t i l /uio and (l 2 Ci)m ax oc Q 2 L 2 , but larger values of L imply 
L > cH~ 1 (z s )/±v t f l /vio and the scaling changes to (l 2 Ci) m ax oc Q 2 . The dotted line (L ~ 80 kpc) shows how power is suppressed in 
scales below the typical bubble size. Panel (b): dependence of the signal on redshift: z = 6.5, 7, 7.5 from top to bottom, corresponding 
to observing frequencies v = 633,593 and 558 GHz, respectively. Panel (c): dependence on the minimum halo mass: M m i n = 5 X 10 5 
(solid line), 5 X 10 8 (dashed line), and 7.5 X 10 8 M Q (dotted line). The distortion y becomes saturated for M min >5 X 10 8 Mq. For 
5 X 10 5 Mq < M m i„ < 5 X 10 8 M , we have the approximate scaling (l 2 Ci) m ax oc Q 2 that is slightly distorted by the presence of the bias 
factor enhancing the dashed curve over the solid one. For M m i n = 7.5 X 10 8 Mq we are well in the optically thick regime (C; oc (7VX 3 ) 2 ) 
and the typical bubble sizes have shrinked dramatically (see Fig. [2] this translates into a clear drop in the power (thinnest curve). 



the redshift: for the range of z s = 6.5, 7, 7.5 (from top to 
bottom) , the amplitude of the power spectrum varies by 
nearly two orders of magnitude. This steep dependence 
arises because the efficiency of the pumping is exponen- 
tially dependent on the difference between the CMB and 
the spin temperatures, and this difference is reduced at 
higher z (see Paper I for more discussion). 

Panel (c) in Figure |4] shows the dependence on the 
minimum size of the halo producing metal-bubbles, with 



M min = 5 x 10 5 (solid line), 5 x 10 8 (dashed line), and 
7.5 x 10 s M Q (dotted line). This roughly brackets the 
possible range of minimum halo sizes, bas ed on the pres- 
ence or absence of H2 molecules (e.g. lHaiman et al.l 
I1997H . At M mm = 5 x 10 5 M , L - cH- 1 {z,)Av th /v 10 , 
which would yield the scaling (l 2 Ci) max oc Q 2 , i.e., inde- 
pendent on M m i n . However, the actual amplitude does 
vary with mass in the range of 5 x 10 5 M© < M m i n < 
7.5 x 10 s M due to several reasons. First, increasing the 
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halo mass results in a larger bias factor (from b ~ 1.4 to 
b ~ 2.6). Second, at large M m i n , there are fewer bubbles 
per unit volume, assumed to contain the same amount of 
metals, so the bubbles become optically thick, and y ~ C? 
and (l 2 Ci) max oc £7 2 (NL 3 ) 2 . But most importantly, for 
M rai „ > 5 x 10 8 M Q , the typical bubble radii (and corre- 
sponding volumes) shrink considerably (c.f.Fig. and 
due to the scaling C; oc (AL 3 ) 2 , this causes a corre- 
spondingly steep drop in the power amplitude (dotted 
line in Fig.(f5J:)). We also note that for the largest mini- 
mum mass considered, the shot noise contribution starts 
dominating at the smallest scales. We find that the am- 
plitude of the signal depends strongly on the parameters 
of the model. However, the toy model for the metal dis- 
tribution was introduced only to gain insight on the na- 
ture of the signal and should be regarded as a toy model. 
Numerical simulations with radiative transfer should be 
used to make more quantitative predictions and will be 
crucial in the interpretation of future observations. 

The multipole corresponding to (l 2 Ci) max in Figure 2] 
is related to the projected k mode at which the convolved 
matter power spectrum peaks. In principle, the typical 
bubble clustering length could also be detected in the 
radial direction if the observing instrument had enough 
spectral resolution: in such case, observations at slightly 
different frequencies would give rise to correlated maps, 
since they would be probing overlapping shells centered 
at similar redshifts. The combination of maps obtained 
at different frequencies could also improve the S/N ratio 
of the final detection. 

Because at the frequencies of interest the signal will be 
dominated by infrared emission from dusty galaxies, it is 
important to be able to remove them from the observed 
map. One way this can be achieved is by obtaining high 
angular resolution images of the sky at the frequencies of 
interest, so that Olber's paradox is avoided and individ- 
ual galaxies can be clipped from the maps. The Atacama 
Large Millimeter Array (ALMA) may be an excellent in- 
strument for this purpose because of its high sensitivity 
and angular resolution (of the order of 0.5"). The draw- 
back is that the field of view of ALMA is small ( Ss 1 ar- 
cmin). However, since our signal peaks at scales of ~ tens 

6 CCAT URL site: |http : //www . submm . org] 



of arcseconds and decays slowly at smaller scales, this 
small field of view should not a fundamental limitation. 
The signal could also be within reach of the sensitivity 
of for thcoming detectors (e.g., SCUBA2, (jHolland et al.l 
120061 )) and of planned single dish experiments such as 
CCAT 6 . However, current small scale CMB experiments, 
with a typical sensitivity of ~ /iK at arcmin scales, are 
still rather far from imposing interesting constraints on 
this effect. 

7. CONCLUSIONS 

In this paper, we have computed the angular fluctu- 
ation of the effect we presented in an earlier work (Pa- 
per I): the CMB spectral distortion induced by 01 63.2 
(itm fine transition, pumped at the epoch of reionization 
by the UV background radiation via the 01 Balmer a 
line. Since the oxygen distribution in space is associ- 
ated with the most over-dense regions hosting the first 
star forming activity at the end of the Dark Ages, there 
must be fluctuations in the angular pattern of the OI- 
induced y distortion. We have considered the cluster- 
ing properties of the signal and computed its angular 
power spectrum using a toy model for the metal distribu- 
tion. We find that the predicted signal is small, but that 
for certain enrichment models it could be detected using 
ultra-deep observations with a sensitive instrument such 
as ALMA, SCUBA2 and CCAT. Even in the event of 
a non-detection, future infrared observations could place 
interesting limits on the metal enrichment of the universe 
during the Dark Ages. This could open the possibility of 
measuring the metal enrichment of the universe before 
it was re-ionized. 01 observations would also be comple- 
mentary to HI 21cm measurements that will map out in 
the near future the distribution of neutral H during the 
epoch or re-ionization. 
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